Ultrathin oxides: bulk-oxide-like model surfaces or unique films? 
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To better understand the electronic and ciiemical properties of wide-gap oxide surfaces at the 
atomic scale, experimental work has focused on epitaxial films on metal substrates. Recent findings 
show that these films are considerably thinner than previously thought. This raises doubts about 
the transferability of the results to surface properties of thicker films and bulk crystals. By means of 
density-functional theory and approximate GW corrections for the electronic spectra we demonstrate 
for three characteristic wide-gap oxides (silica, alumina, and hafnia) the influence of the substrate 
and highlight critical diff'erences between the ultrathin films and surfaces of bulk materials. Our 
results imply that monolayer-thin oxide films have rather unique properties. 

PACS numbers: 68.47.Gh; 73.20.-r; 71.20.Ps; 

Keywords: thin oxide films, electronic structure, silica, alumina, hafnia 



o 



o 

O 

o 



X 



On the nanoscale materials often reveal new and unex- 
pected features. In particular ultrathin oxide films have 
recently generated increasing interest. First, thin oxide 
films themselves play an important role in technologi- 
cal applications such as electronic devices, fuel cells, gas 
sensors, corrosion and scratch protection, or heteroge- 
neous catalysts. Second, they are increasingly used as 
model oxides in the surface science approach to thicker 
or bulk-like oxides: Many experimental techniques such 
as scanning tunneling microscopy or photocmission spec- 
troscopy cannot be directly applied to insulating mate- 
rials, but require electrically conducting samples. Ultra- 
thin epitaxial oxide films on metal substrates overcome 
this limitation [H, 0, i, i, H @, 0, i, i, 113, IH [H • How- 
ever, the direct transfer of experimental conclusions from 
these films to bulk oxide surfaces relies on the two criti- 
cal assumptions that a) the thickness dependence of the 
film and b) substrate-film interactions are negligible for 
the surface properties (cf. Fig. [1]). In the present let- 
ter we argue that the finite thickness and the presence 
of the substrate modify the film's geometry and/or sto- 
ichiometry in a way that is unknown and maybe even 
impossible for surfaces of thicker films and bulk crys- 
tals. In addition, the coupling of the film to the metal's 
Fermi level may induce charge transfer 0. The result- 
ing mechanical, electronic, and chemical properties may 
then differ considerably from the systems that were the 
original goal of the study. Instead, these ultrathin films 
should be considered as new materials with interesting 
and novel properties of their own. 

For a meaningful comparison between thin films and 
bulk oxides one has to first identify the oxide surfaces 
with the closest resemblance to the ultrathin film. Then, 
one would like to isolate the influence of the substrate 
from that of the finite thickness or the thickness depen- 
dence in the surface science approach, to analyze each 
contribution separately. While the latter is not possi- 
ble with the experimental tools available to date, it is 
easily achieved in theory by the hypothetical two-step 
pathway shown in Fig. [TJ a) reduce the thickness of free- 
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FIG. 1: (Color online) The critical assumptions in the surface 
science approach (here for Si02/Mo) are tested by following 
a theoretical pathway that transforms the surface of a bulk 
material (top right) to an epitaxial thin film on Mo(112) (top 
left): a) decreasing the thickness and b) placing the film onto 
the substrate. O atoms are depicted in red. Si in yellow and 
Mo in gray. 



standing films and study the effect of the thickness vari- 
ation b) place an ultrathin film on a substrate analyz- 
ing the substrate dependence. In this Letter, we employ 
density-functional th eory (DFT) in the local-density ap- 
proximation (LDA) [iJI to address these questions for 
three representative thin-film oxide systems: 1) silica on 
Mo(112), a recent prototypical surface science model for 
ceramic support materials of heterogeneous catalysts, 2) 
alumina, which in various modifications is also frequently 
applied as model substrate for catalytic reactions, and 3) 
hafnia, the material of choice in the most advanced mi- 
croprocessor technology . 

For silica, flat and well-ordered films have successfully 
been grown on Mo(112) by Freund et al. and by Good- 
man et al. [3, 07 01 ■ Recent ah initio simulations [3, 0| 
suggest a two-dimensional network structure (isostruc- 
tural, as we will show below, to the a-quartz (0001) 
surface (lit), whereas Goodman et al. favor an iso- 
lated [Si04] cluster model 0]. Indications for thickness 
effects were given by Wendt et al., who have demon- 
strated in a comparative study including thicker (but 
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amorphous) films that the electronic structure varies with 
the film thickness Alumina films have been grown 
on a wide variety of substrates, e.g.. Re. Mo, NiaAl, 
NiAl, or FC3AI. A typical, intensively studied example 
is NiAl(llO), where a well-ordered alumina film of de- 
fined thickness is obtained by a direct oxidation of the 
alloy's surface • Its complex structure is governed by 
the interplay between the local, preferentially hexagonal 
ordering of the oxide and its linkage to the Al atoms at 
the interface to the metal alloy [lH . Recent STM experi- 
ments for the adsorption of small gold cluster on this film 
have shown that surface properties are modified by the 
film's substrate Ultrathin hafnia films on metallic 

substrates have not been investigated so far. 

To investigate how far ultrathin films could refiect 
the surface properties of macroscopic samples in prin- 
ciple we first exclude the substrate influence by study- 
ing free-standing silica, alumina, and hafnia films. For 
all three oxides we find substantial differences between 
monolayer (ML)-thin and thicker films. We concentrate 
on the atomic structure and the electronic density of 
states (DOS). Many other properties are correlated to 
these two key aspects, e.g., the atomic configuration and 
the local electronic structure determine the chemical ac- 
tivity. The films are fully relaxed within the surface unit 
cell of the most stable bulk surface, i.e., the (0001) sur- 
faces of a-quartz (silica) [11] and a-alumina [l3|, and 
the (111) surface of cubic hafnia This ensures that, 
upon increasing the thickness, the films systematically 
develop towards the bulk limit. 

For silica, we have compared a variety of surface ter- 
minations and reconstructions. We find that the most 
stable films are stoichiometric and have fully saturated 
bonds. The structures represent perfect networks of 
[8104] units, in agreement with the trends in bulk sili- 
cates [1^. This similarity in the local geometry is con- 
trasted by the flexibility of the bonding network that 
easily adjusts to the constraints imposed by the finite 
thickness and leads to quite different structures for each 
thickness. Films with more than 4 formula units are best 
described as a-quartz-like layers sandwiched between two 
reconstructed surfaces with three-membered silica rings 
at the interface (SisOio being highly distorted). In 
the Si408 film, which has no intermediate layers, the re- 
constructed surfaces are directly linked together. Si2 04 
consists of a single-layer network with two-membered sil- 
ica rings, a motif known from other silica surfaces (7l.[20|. 

The structural variations in the alumina and hafnia 
films are less drastic, but significant for ultrathin films, 
too. For example, the outermost Al atom of the alumina 
surface is known to strongly relax inward 17|. This re- 
laxation shows a pronounced thickness dependence. The 
interlayer spacing is reduced by 62% for the AI4O6 film, 
by 96% for AlgOg, and 87% for AlgOia and thicker films. 
The relaxations at the hafnia surface are smaller (<4% of 
the interlayer spacing), but the thinnest films again de- 
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FIG. 2: Surface DOS for silica, alumina and hafnia films in 
comparison to the bulk DOS. 



viate from their thicker counterparts. The trend towards 
smaller structural changes reflects the higher coordina- 
tion in alumina (6/4 for Al/0, respectively) and hafnia 
(8/4 for Hf/0) compared to sifica (4/2 for Si/0) which 
limits the geometrical flexibility. 

The effect of the finite thickness on the electronic and 
chemical properties can be analyzed in terms of the DOS, 
displayed for selected films in Fig. [2] While the wide-gap 
insulating character of the perfect films is not affected 
by the thickness, the peak structure of the DOS differs 
significantly for ultrathin films. Since defect levels are 
most likely subject to the same shifts as the bands from 
which they derive, this implies that ML-thin films may 
show critical deviations in their defect properties. Fig. [2] 
highlights the differences in the local DOS at the surface 
(the detailed evolution of the electronic structure will be 
discussed elsewhere) comparing the surface DOS of one 
or two ultrathin films to a thick film, more representative 
for macroscopic samples, and the total DOS of the bulk. 
In all three cases, the ultrathin films have a unique peak 
structure that differs from that of thicker films. These 



3 



varying peaks are not associated vv^ith particular local 
orbitals but must be considered properties of the overall 
electronic structure. In addition, the AI4O6 film shows a 
remarkable reduction of the band gap, while the valence 
band width is comparable to that of thicker films and 
also the bulk. This indicates an overall weakening of the 
Madelung potential when a substantial amount of the 
ions is undercoordinated [l3|. In all other cases, the band 
widths and gaps are close to their bulk values. In general, 
we find that the changes in the electronic structure follow 
those in the local atomic structure, which underlines the 
crucial importance of the latter. 

Of course, supported films are additionally infiuenced 
by the interface to their substrate. Yet, our results for 
free-standing films clearly demonstrate that the ML-thin 
films cannot be considered representative for the struc- 
tural, electronic, defect, or surface properties of the re- 
spective bulk materials in contrast to earlier statements 
[J, 0, [3, [l^l. Only when the films are thicker than 3-4 
layers, the surface properties of defect-free surfaces are 
approaching those of bulk-like systems. This is a no- 
ticeable smaller film thickness than known for metal sur- 
faces, where for closed packed surfaces typically around 
8 in some cases even up to 20 layers [ilj are needed, but 
in line with other ionic materials such as NaCl [i^]. We 
attribute this trend to the degree of localization of the 
valence electrons in these materials, which is larger in 
the insulators (considered here) than in metals with typ- 
ical semiconductors in between. The deviating behavior 
of ML-thin films offers the prospect that experimentally 
accessible ultrathin oxide films such as silica on Mo(112) 
may be understood as new and unique materials in their 
own right, opening possible new routes to devise insula- 
tor surfaces with novel properties. 

We will now in more detail discuss the role of the sub- 
strate for the example of silica on Mo(112), assuming 
the recently suggested structure model 0, 01 (denoted as 
"siloxane surface"): Corner-sharing [Si04] tetrahedrons 
form a two-dimensional, perfectly bonded hexagonal net- 
work, that is linked to the Mo surface via oxygen (cf. 
Fig. [T]). The film has a Si205 stoichiometry that does 
not naturally occur in bulk silica. A comparison to other 
known thin film and surface structures reveals that it is 
identical to the "dense reconstruction" of the a-quartz 
(0001) surface [H], the "silicate adlayer" on 4H and 6H 
SiC [2^ , and the "siloxane surface" of clay minerals . 
The correspondence to the a-quartz surface (which ap- 
parently has not been recognized previously) makes the 
silica film on Mo(112) an ideal candidate for studying the 
role of the metal substrate. For this, we regard the a- 
quartz layers below the siloxane surface as an alternative, 
insulating substrate in place of Mo(112). 

Focusing on the structural aspects first, we find that 
the siloxane surface adapts to the substrate by reorient- 
ing the [8104] tetrahedrons as rigid units, a behavior 
well-known for silicates [31 . A substrate-independent. 
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FIG. 3: Substrate influence on the silica surface DOS for a) 
the Mo(112) and b) the Q-quartz substrates in comparison 
to the total DOS of bulk a-quartz c). The bottom of the 
surface conduction band is marked by an arrow. The energetic 
position of the Mo d-bands is indicated by the gray bar. 



"ideal" siloxane surface structure was determined from 
a free-standing siloxane double layer (Si408). Its opti- 
mized lattice constant (5.24 A) is very close to the quasi- 
hexagonal c(2 X 2) Mo(112) superlatticc (a=5.2lA and 
a=63°) and also to the (V3xv^)R30° pattern (5.25 A) 
observed for SiC [2^, whereas the lattice mismatch to 
a-quartz (4.85 A) amounts to 8%. This provides a sim- 
ple explanation why the siloxane surface forms readily on 
the lattice- matched substrates Mo(112) and 4II/6H-SiC 
(0001), while it has not been possible to prepare it ex- 
perimentally on a-quartz. A further important factor is 
the linkage of the silica film to the Mo rows on Mo(112). 
The threefold symmetry axis of the [Si04] tetrahedrons is 
tilted by 14° against the surface normal (cf. Fig.[T]). This 
introduces a height difference (buckling) for the surface 
atoms of ~0.4A. The siloxane surface on quartz, how- 
ever, tilts only by 4° and buckles by 0.2 A. We conclude 
that the siloxane surface on Mo(112), used in the surface 
science experiments, exhibits new features and deviates 
noticeable from the ideal a-quartz surface despite the 
identical chemical connectivity. 

To visualize the substrate's influence on the electronic 
structure, the local DOS for the siloxane surface on 
Mo(112) and a-quartz are presented in Fig. [3] along with 
the total DOS of bulk a-quartz as a reference. The 
band edges at the surface agree with those of a-quartz to 
within 0.5 eV for both substrates. There are no surface 
states in the large band gap. The Mo metal states, which 
energetically lie within the surface band gap, play no im- 
portant role at the surface. However, the peak structure 
in the valence region of the surface DOS differs for the 
two substrates and deviates also from the DOS of bulk 
a-quartz. We have analyzed the origin of these changes 
by placing the silica/Mo surface structure on quartz and 
find that mainly the structural variations are responsible 
for the changes observed, whereas the chemically differ- 
ent interface plays a minor role. 

That the surface electronic structure of the well- 
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TABLE I: Peak positions relative to the Fermi level for the 
electronic surface DOS of the siloxane surface on Mo(112) 
in comparison with experiment. The energy scale is shifted 
by 5.1 eV compared to Fig. [3] The theoretical data are not 
corrected for the estimated GW shift of 2eV (see text). 



pcrimentally accessible silica film on Mo(112). We expect 
that similar conclusions hold for other metal-supported 
oxide films. Thus, ML-epitaxial thin films should not be 
seen merely as models for thick oxides in heterogeneous 
catalysis or other applications. Instead they may offer a 
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ordered silica film on Mo(112) possesses a unique peak 
structure has been observed in ultraviolet photoelec- 
tron spectroscopy (UPS) and metastable impact electron 
spectroscopy (MIES) 0, 0. For comparing the DFT- 
LDA DOS to experimental spectra, we focus on the peak 
positions since the spectral intensities depend on the 
spectroscopic method and even make some peaks dis- 
appear in the experimental MIE and UP spectra 
To overcome limitations of DFT-KS in describing the 
quasiparticle spectrum probed in experiment, we employ 
many-body perturbation theory in the GW approxima- 
tion 25|, l26[. In GW, the electronic self-energy (which 
connects the non-interacting electrons of DFT-KS to the 
truly interacting ones) is given by the product of the 
Green function G and the screened interaction W. GW 
is typically applied as a correction scheme to the DFT- 
LDA energies which corrects deficiencies of the LDA and 
includes quasiparticle effects 25,2^. We have performed 
GW calculations for bulk a-quartz. Transferring in a 
first crude approximation the bulk corrections to the sil- 
ica film, we estimate that the silica features in the valence 
region are shifted downward by ~2 eV. Taking this GW 
correction into account, the agreement with experiment 
(cf. Tab. HI becomes very good. This agreement not only 
supports the siloxane surface as structural model for the 
silica film on Mo(112), but also shows that the peaks ob- 
served in experiment arc not generic silica features, but 
specific to the actual film. 

In light of these results we argue that the silica film 
on Mo(112) must in fact be considered a novel oxidic 
surface, as envisioned from free-standing films, and is not 
simply another incarnation of a siloxane surface on bulk 
a-quartz. The lattice mismatch of 8% to a-quartz (which 
may prevent the formation of large ordered domains in 
real samples) induces important changes in the structure 
compared to Mo(112), which subsequently modifies the 
electronic structure. 

In summary we find that ultrathin films of wide-gap 
oxides require two or three bulk-like layers to develop the 
characteristic properties of the bulk materials and their 
surfaces. This implies that ML-films are no simple ana- 
log of thicker films or realistic oxide surfaces. In addition, 
we demonstrated significant substrate- induced changes in 
the surface structure and electronic properties for the ex- 
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